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A study of magneto-transport through quantum dots is presented. The model allows to analyze 
tunnelling both from bulk-like contacts and from 2D accumulation layers. The fine features in the 
I-V characteristics due to the quantum dot states are known to be shifted to different voltages 
depending upon the value of the magnetic field. While this effect is also well reproduced by our 
calculations, in this work we concentrate on the amplitude of each current resonance as a function 
of the magnetic field. Such amplitudes show oscillations reflecting the variation of the density of 
states at the Fermi energy in the emitter. Furthermore the amplitude increases as a function of 
the magnetic field for certain features while it decreases for others. In particular we demonstrate 
that the behaviour of the amplitude of the current resonances is linked to the value of the angular 
momentum of each dot level through which tunnelling occurs. We show that a selection rule on the 
angular momentum must be satisfied. As a consequence, tunnelling through specific dot states is 
strongly suppresses and sometimes prohibited altogether by the presence of the magnetic field. This 
will allow to extract from the experimental curves detailed information on the nature of the quantum 
dot wavefunctions involved in the electronic transport. Furthermore, when tunnelling occurs from 
a 2D accumulation layer to the quantum dot, the presence of a magnetic field hugely increases the 
strength of some resonant features. This effect is predicted by our model and, to the best of our 
knowledge, has never been observed. 



I. INTRODUCTION 

The manufacturing of structures where the electronic 
motion is confined in 1,2,3 directions has been possible 
for several years and a large amount of experimental 
evidence has been collected to study the behaviour of 
electrons in 2,1,0 dimensions. This has resulted in a se- 
ries of new physical phenomena and has raised hopes of 
building devices where the behaviour of single electrons 
is controlled. Quantum dots are structures where the 
electron motion is confined in all directions and there- 
fore their energy spectrum is discrete. Several spectro- 
scopic techniques such as linear and non-linear trans- 
port, capacitance measurements and far-infrared spec- 
troscopy have been used to investigate the energy spec- 
tra of quantum dots. The addition of an external mag- 
netic field has proved to be an excellent way to probe the 
low-dimensional systems since it introduces a structure- 
independent quantisation that is superimposed to that 
induced by the fabrication process. Quantum dots ob- 
tained from double-barrier resonant tunnelling semicon- 
ductor heterostructures are now used by several groups 
for magneto-tunnelling investigations [|l]-Q, since it is 
possible to probe the quantum dot spectra in several 
regimes from the very low to the very high electron num- 
ber occupancy, tuning the relative strength between con- 
finement and electron-electron interactions. The spec- 
troscopic investigation is still exciting and far from com- 



plete, and several issues such as the importance of 
the electrostatic interaction among the carriers, detailed 
characteristics of the quantum dot spectra, degeneracy of 
the states and spin of the different multiplets must still 
be fully solved. 

Practically, in all the previous studies on non-linear 
transport, the emphasis was placed on the displacement 
of the I-V features produced in the quantum dot states 
as a function of the magnetic field in order to probe the 
quantum dot spectrum On a more fundamental 

point lies the issue of actually probing the confined wave- 
functions inside a quantum dot. This analysis has been 
carried out by P.H. Beton et al. || for the states in a 
quantum wire. The method consists in measuring the 
current voltage characteristics of a quantum wire in the 
presence of a magnetic field as in the previous studies, 
but to focus the analysis on the magnitude of the cur- 
rent of each peak. Indeed, it is shown that the amplitude 
of the current can be associated with the Fourier trans- 
form of the electron wave function in the direction across 
the wire. Very recently we reported on a similar 
analysis of the experimental data collected on quantum 
dots fabricated from a double-barrier resonant tunnelling 
structure by an ion implantation and electrode regrowth 
process Rj. To the best of our knowledge this was the 
first time such an analysis had been performed for quan- 
tum dots. 

In this work we lay the theoretical foundations to ob- 
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tain a detailed understanding of the electronic proper- 
ties of quantum dots from magneto-transport experi- 
ments. In particular we address the problem of non- 
linear magneto-transport in 3D-0D-3D and 2D-0D-2D 
structures. This can yield informations not only on the 
energy position of the dot states, but also on the ac- 
tual properties of the wavefunctions, i.e their angular 
momentum and their spatial extension. Extending the 
work performed by P.H. Beton et al, we find that also for 
quantum dots the height of the resonances in the current- 
voltage characteristics depends strongly on the magnetic 
field. The height of each resonance shows oscillations 
as a function of the magnetic field which reflect the os- 
cillations in the density of states at the Fermi energy. 
Furthermore the behaviour of the amplitude of the fine 
features changes considerably depending on the particu- 
lar quantum dot state considered. It is found that a gen- 
eral trend is obeyed depending on the value of the angular 
momentum of the dot state, and that the angular momen- 
tum introduces a selection rule that strongly suppresses 
tunnelling through states with angular momentum co- 
linear to the magnetic field. Above certain values of the 
magnetic field the selection rules totally prohibits tun- 
nelling through such states. The last section deals with 
the 2D-0D magneto-transport occurring in those experi- 
mental structures where an accumulation layer is formed 
in the emitter. Here the effect of the magnetic field is 
totally different from the 3D-0D-3D case. We show in 
particular that the amplitude of the current resonances 
are a direct measure of the overlap between the Landau 
levels in the emitter and the dot states involved in the 
tunnelling process. The presence of the magnetic field 
hugely increases the tunnelling current flowing through 
the quantum dot ground state and indeed through all 
the states with principal quantum number n=0 and non- 
positive angular momentum. It must be noted that this 
effect applies to the quantum dot ground state, i.e. when 
only a single electron occupies the dot. Therefore it re- 
mains valid even in the presence of electron-electron re- 
pulsion inside the dot. The theoretical analysis developed 
here will help to distinguish the nature of the various fea- 
tures in the experimental investigations ||, to establish 
unambiguously the dimensionality of the electrodes and 
to gain further knowledge of the properties of quantum 
dots. 

II. GENERAL REMARKS 

The current-voltage characteristics of quantum dot 
structures are dominated by the concept of resonant tun- 
nelling. The current flows through a specific quantum dot 
state only if its energy falls within the energy range of 
the occupied states in the emitter conduction band. The 
total current flowing through the structure for a given 
applied voltage bias is the sum of the current flowing 



through each quantum dot state. In systems such as 
double-barrier resonant tunnelling heterostructures, at 
zero voltage no current flows through the structure be- 
cause the quantum dot states lie above the Fermi energy 
in the emitter. As the voltage is raised, the quantum dot 
states fall below the Fermi energy, resonant tunnelling 
becomes possible and the current increases. It follows 
that the voltage at which each dot state becomes avail- 
able for tunnelling is associated with the energy of the 
dot state itself. Therefore transport measurements pro- 
vide informations on the quantum dot spectrum. 

It has long been argued that transport measurements 
must be performed in the presence of a magnetic field in 
order to be a powerful spectroscopic tool. This is because 
the presence of impurity states in the quantum dot region 
can produce features in the I-V characteristics that are 
similar to the ones produced by the quantum dot states. 
The presence of a magnetic field allows to alter the geom- 
etry and the strength of the confinement of the quantum 
dot in an controllable fashion so that the properties of 
the quantum dot can be differentiated from those due 
to impurities. Furthermore the features produced by the 
quantum dot can be analyzed in much greater details 
by the application of a magnetic field, hence the impor- 
tance of understanding in details the physics of magneto- 
tunnelling in quantum dots. 

The typical structure to study vertical transport in 
quantum dots consists of a double-barrier quantum well 
which is confined in the lateral directions. The lateral 
confinement can be introduced by means of several meth- 
ods such as etching QUII , ion-bombardment [j7| , and the 
presence of gates pj . For all these methods the resulting 
confining potential can be very accurately approximated 
by a parabola. In the growth direction the heterostruc- 
ture confining potential V(z) is much stronger than the 
lateral one and has the shape of a square well. Therefore 
the Hamiltonian for the dot containing N particles can 
be written as: 

N ( - X \2 N 

H = E - 2 m + 1™*" 2 * S>? + yf) + n*) a) 

i=l i=l 

In this work we are going to analyze the case where 
the magnetic field is applied along the vertical direction 
of the structure, i.e. along the direction of the current. 
As we are going to deal with small structures where the 
single electron effects dominate over the electron-electron 
repulsion ||, we do not take into account the Coulomb 
repulsion inside the dot. This approximation consider- 
ably simplifies the problem of transport through quan- 
tum dots in the non-linear regime and, while it is not the- 
oretically satisfactory, it has proven correct in several cir- 
cumstances both experimentally |l| [||]]6]|| an d theoret- 
ically [p| ]Kl JTT| ] . Nevertheless, the analysis presented here 
is valid also in the presence of the electron electron inter- 
action when this can be approximated with the Hubbard- 



2 



like hamiltonian. This has been proved a valid approach 
for small dots of radius 50 nm or less || , as those ones 
considered here. While many theoretical efforts, based on 
non-equilibrium Green's functions, have been devoted to 
the development of a theory of transport in the non- linear 
regime for correlated systems ]l2| |, a comprehensive the- 
ory capable of dealing with a multi-electron system and 
including partial occupancy and non-equilibrium filling 
of the dot states in a quantitative fashion does not exist. 
In the single particle approximation eq.(l) can be solved 
exactly. Thus the single electron spectrum, neglecting 
the effects of spin, is given by: 

E nd ,m d = hu d (2n d + \m d \ + 1) + ^-fruj c m d + E (2) 



where ui d = y uJq + , u>q characterizes the strength 
of the lateral confinement in the dot, u c is the cyclotron 
frequency and Eg is the energy of the level due to the con- 
finement in the growth direction. The principal quantum 
number n d can take any non-negative integer value. On 
the contrary, the quantum number m d , associated with 
the angular momentum in the direction of the current 
and the magnetic field, can take any positive or negative 
integer value. The eigenfunctions related to the eq. (2) 
are given by: 

$ d (r, ip, z) = -—R nd , md {r)* d {z) (3) 

V 27T 

where R nd , md (?*) is a function that can be expressed in 
terms of the hypergeometric function E^ , E3| : 

Rn, m (r) oc e-« /2 e |m|/2 F(-n, |m| + 1, £) (4) 

where £ = uj d ra*r 2 /h. ^ d (z) is the wavefunction in the 
growth direction which consists of the familiar ground 
state of a square well of finite height. 

As pioneered by Bardeen [|l5| and further developed 
by Payne jl6) and by Liu and Aers the tunnelling 
current can be calculated in the framework of a trans- 
fer matrix approach. The system is divided into three 
sub-systems, the two contacts and the quantum dot, and 
the spectrum for each sub-system considered isolated is 
determined. 

The tunnelling process is modeled as a two-step pro- 
cess in which electrons first tunnel from the emitter to the 
dot and eventually from the dot to the collector. In spite 
of its sequential nature, this approach is well-known to 
reproduce the correct physics even in the resonant tun- 
nelling regime pqj. The current for each step can be 
calculated from the knowledge of the wavefunctions in 
the three different regions. At equilibrium the current 
flowing from the emitter to the dot must be equal to the 
one going from the dot to the collector. Imposing this 
condition, enables one to calculate the current as well 
as the occupancy of the dot states. The strength of the 



tunnelling current flowing through each dot state is deter- 
mined by the overlap integral between the dot wavefunc- 
tions and those of the contacts ]l^ , |l7| ,p|,^9| . Following 
Bardeen, this is given by: 

Me^d = (J^j ^($eV$5 - n^ e ).dS (5) 

for the transition of the electrons from the emitter to the 
dot, here S is any surface across the emitter-dot barrier. 
A similar formula applies for the tunnelling between the 
dot and the collector. 

In order to calculate the current it is necessary to know 
not only the electron wave functions in the dot, but also 
those of the contacts. As these depend on the actual 
shape of the contacts, we are going to discuss two differ- 
ent cases separately. 

III. FIRST CASE: 3D CONTACTS 

Let us first suppose that no geometrical confinement is 
present in the contacts. This is indeed the case in many 
experimental situations f|,f|;|]j7|] where dots are either ob- 
tain by means of deep etching or by ion bombardment. 
In the presence of an applied magnetic field, the energy 
spectrum in the contacts takes the form of highly degen- 
erate ID Landau sub-bands which are occupied up to the 
Fermi energy. The energy spectrum is given by: 

^ *. / m e + \m e \ 1. 7i 2 fc 2 

E ne , me ,k z = Ku c (n e + 2 ' e| + -) + ^-i (6) 

where k z is the momentum in the growth direction. As 
for the quantum dot, m e can take any integer value while 
n e can take all non-negative integer values. Note how- 
ever the difference between positive and negative values 
of m e . While each Landau sub-bands contains all nega- 
tive values of m e , the lowest sub-band, of minimum en- 
ergy ^hijj c , does not contain positive values of m e , the 
first excited sub-band contains only m e — 1, the second 
excited sub-band m e = 1, 2 and so on.. In general terms 
the maximum value of the angular momentum in the j-th 
sub-band is m e = j — 1. 

The eigenfunctions are plane waves along the vertical 
growth direction while for the lateral part they are of the 
same shape as those of the quantum dot: 

<f> e (r,tp,z) = — =R ncime (r)^ e {z) (7) 
V 2n 

where ^ e {z) is the wavefunction defined by the one di- 
mensional square potential: V(z) — before the emitter 
barrier for z < b and V(z) = Vq for z > b, b being the 
abscissa of the left side of the barrier and Vq the height 
of the barrier. In this case ^ e (z) is given by: 
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n(k e L e )evp(-K e [z - b}) 



z<b 
z>b 



(8) 



where L e is a normalization constant, k e = \/2m*E/fr 
is the lateral momentum and K e is defined as K e — 
y/2m*(Vo — E)/fi. E is the energy in the vertical direc- 
tion. The spectrum and wavefunctions of the collector 
are exactly of the same shape as those of the emitter. 

The overlap between the quantum dot states and those 
in the contacts depends on the height and width of the 
barriers along the direction through which the current 
flows. In the lateral directions it depends on the particu- 
lar quantum dot and contact states considered. The over- 
lap integral for a one dimensional problem with barriers 
of square shape has been studied in details by Payne jl6) . 
In this work we will adopt the same approach to treat 
the overlap in the growth direction and we will extend 
the method to the three dimensional case. In order to 
do so one must calculate the overlap integrals in the lat- 
eral directions between the quantum dot and the contact 
states. These are of the form: 
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where the e and d refer to states in the contacts and in the 
dot respectively. The overlap integral || can be calculated 
analytically for all values of rrid, rid, m e , n e . Furthermore 
it is clear from ^| that only the states with rrid = rn e = m 
produce a non zero overlap. Therefore from the highly 
degenerate Landau sub-bands in the contacts, at most 
only one particular state will contribute to the tunnelling 
current flowing through each dot state: the one with the 
right value of the angular momentum rrid- For example 
from the lowest emitter sub-band there cannot be any 
contribution to the tunnelling through those dot states 
with positive angular momenta. 

The overlap integrals || can be calculated analytically 
and are given by [ jl3| : 



s mi , me ^r( 7 )A- n *- n «-i f (A - k d y«(A - k e y 



XF(-Ud, -»<s,7, 



kdk e 



(\-k d )(\-k e ) 



) 



where 



&d,e 



(|"d,e|+"d,e)! 



ad,, 



2m'LO diC {B) ■ 



X 



1 <+< U J _ 



2 |m <J.^n d , e !|m d , e |! 2 

F is the hypergeometric function [ p^[ . 

7 = \md, e \ + 1 and uj e = uj c j2. 



4 a 2 d xa 2 e ' , " a < e ~ 2a* c 

Therefore combining formula ^ with the overlap in- 
tegrals M in the direction of the current as outlined 
by Payne and by Liu and Aers [^7), one can calcu- 
late the current flowing through the quantum dot in an 



analytical way. The above technique has been adopted 
to calculate the I-V characteristics of 3D-0D-3D struc- 
tures in the presence of an external magnetic field ap- 
plied in the vertical direction, i.e. parallel to the current. 
Fig. [I] shows the calculated I-V characteristics at T = AK 
for several values of the magnetic field, calculated with 
Hluq = 15meV. This value of the confinement strength 
together with our choice of Gao.33Alo.Q7As barriers of 
widths 8.7nm and of a GaAs quantum well 5.1nm width, 
is in accordance with some of our experimental devices 
and measurements At low temperature as the volt- 
age bias is raised, electrons tunnel towards unoccupied 
states in the collector so that tunnelling is dominated by 
the emitter barrier. 

The triangular shape of the peaks at B=0T is typical 
of a 3D-0D tunnelling. As expected, the presence of the 
magnetic field shifts the position of the peaks on the volt- 
age scale. This effect is due to the change in the energy of 
the quantum dot states as the magnetic field is increased. 
More important is the current amplitude of the peaks 
that oscillates as a function of B and have a different be- 
haviour from peak to peak. To better illustrate this effect 
in fig. ^ we show the behaviour of the current amplitude 
for different states. States rid = 0, = —2,-1,0,1,2 
and n<j = I, vrid = are represented. 

The amplitude for all the states shows oscillations as 
a function of the magnetic field. These oscillations have 
the characteristic dependence and are due to the os- 
cillations of the density of states at the Fermi energy in 
the emitter as the magnetic field is changed. 

Furthermore, as the magnetic field is increased, the 
peaks show a rather different behaviour in terms of their 
amplitude. The peak associated with the ground state 
shows an amplitude that decreases systematically as a 
function of the magnetic field, while the first excited state 
rid = 0, rrid = — 1 initially increases in amplitude, attains 
a maximum, and eventually decreases. The state rid — 
0, rrid — —2 has the same qualitative features as the rid — 
0, rrid = — 1 peak, but the maximum is more pronounced 
and shifted at higher magnetic field. These effects have 
indeed been observed recently by B. Jouault et al 

The presence of the magnetic field produces two effects. 
Firstly it reduces the momentum along the direction of 
the current for the electrons in the contacts by moving 
the sub-bands' minima to higher energies. Secondly it 
confines the electrons towards the center of the dot. The 
first effect tends to diminish the current, while the second 
one increases it by increasing the overlap in the lateral 
direction between the quantum dot wavefunctions and 
those of the contacts. Thus each Landau sub-band in 
the emitter gives a current which attains a maximum 
at a finite magnetic field. The position of this maximum 
depends on the relative strength of the two effects. More- 
over increasing the angular momentum nid decreases the 
overlap integral and shifts this maximum to higher fields. 
The total current is obtained by summing the contribu- 
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tion of all the occupied Landau levels in the emitter and 
it follows that the shift of the maximum current peak 
from the origin B — OT is a direct measurement of the 
magnitude of the angular momentum mj. 

Perhaps even more striking is the totally different be- 
haviour of the states rid = 0, rrid = 1 and rid = 0,rrid — 
— 1. At B = OT these states are degenerate and the ra- 
dial part of the wave function is identical for both of them 
for all values of the magnetic field. This is confirmed by 
the fact that at B = the amplitude is the same for 
both of them. However, the lowest Landau sub-band in 
the contacts does not contain the state rid = Q,md = 1 
so that there is no contribution to the current flowing 
through such state from this Landau sub-band. Since 
the radial part of the overlap integrals is at a maximum 
between states with the same principal quantum number, 
as soon as well defined Landau levels are formed in the 
contacts tunnelling through the state rid — 0,rrid = 1 is 
severely limited. Indeed whenever the magnetic field is 
strong enough so that only the lowest Landau sub-band 
is occupied in the contacts, tunnelling through the states 
with angular momentum to > is altogether forbidden. 
This is a manifestation of the general selection rule previ- 
ously mentioned, tunnelling through a state with positive 
angular momentum to is forbidden whenever there are 
only to occupied Landau sub-bands in the contacts. The 
consequence of this is clearly observed in fig. [l]: all the 
current peaks which are not related to the first Landau 
level in the emitter disappear quickly when the magnetic 
field increases. These peaks related to a positive angular 
momentum with n = or to a strictly positive value of 
n (with all values of to) shift to higher voltage bias as 
the magnetic field increases. As an example the second 
current peak due to the first excited quantum box state 
at about 15meV splits into two peaks with B (n = 0, 
to = ±1), the third one into three (n = 0, to = ±2 and 
n = 1, to = 0) and so on. All the current contributions 
of the peaks shifting to higher voltages (n — 0, to > or 
n > 0) wash out very quickly with the magnetic field. 

Finally fig. || shows that the contribution of the first 
Landau level for the state rid = lj^d = is weak rel- 
atively to the other dot states rid = 0,TO(j = 0,-1,-2. 
This effect reflects the poor overlap of the radial part be- 
tween the states n e = 0, m e = and rid — ljTO^ = 0. 
Therefore, the current is sensitive to the shape of the 
confining potential, even if the shift of the maximum is 
a general feature for every type of quantum dot. 



IV. SECOND CASE: THE 2D-0D TRANSPORT 

The same calculations can be performed for a one- 
dimensionally confined emitter, which corresponds to 
those experimental structures where an accumulation 
layer is formed in the emitter. In this case a two- 
dimensional electron gas (2DEG) is formed to the left- 



side of the emitter-dot barrier and the transport proper- 
ties of the structures are dominated by the resonant tun- 
nelling between the 2DEG and the quantum dot. The 
form of the confining potential that produces the 2DEG 
in the contacts is roughly of triangular shape and only the 
lowest level E® of this accumulation layer is usually filled 
and participate to the tunnelling current p(J . Therefore 
we are only going to consider the lowest emitter state E® 
in the following analysis. The method adopted in the pre- 
vious section for the calculation of the I-V characteristics 
can be readily used also in this case. The only difference 
is that the energies and eigenfunctions of the emitter are 
modified and can be written as: 

„ , TO„ + ItoJ 1 , n , 

E ne , me = huj c (n e + li e| + -) + E° e (11) 

the wavefunctions are always given by the equation |^ 
but here ^ e (z) is the same for all the electrons in the 
accumulation layer. In other words the presence of a 
magnetic field produces a fully-discrete spectrum in the 
accumulation layer. Several different forms of the elec- 
tron wavefunction can be taken in the triangular groove. 
However for the scope of our calculations the detail of 
the wavefunction is not important, the only restriction 
being that under the barrier it attains an exponentially 
decaying form. Furthermore as such wavefunction is the 
same for all electrons, the overlap integral in the growth 
direction is a constant. Assuming a simple Lorentzian 
spectral function of characteristic width V for each Lan- 
dau level in the accumulation layer [2lf| , the current can 
be calculated as: 



I = 2e\M 



±\2 



E 



M, 



n d,' m d n e 



f(E d - E f )T/n 



n d, rrL d :n e ; m e 



^Eld 



r 2 



(12) 



where A.E e ^d is the difference in energy between the state 
in the emitter and the state in the dot at a given voltage 
bias. 

Figure || shows the contribution to the I-V characteris- 
tics for the ground state in the dot as the magnetic field 
increases. At B=0T the ground state current increases 
exponentially with the energy, as expected if there is nei- 
ther magnetic nor geometric lateral confinement in the 
contacts. Indeed the overlap integral is greater with the 
states of the emitter near the conduction band edge be- 
cause their small lateral momentum is of the same order 
of magnitude than size of the dot state. Thus the sharp 
closure at V=20meV arises when the bottom of the con- 
duction band in the emitter is aligned with the state of 
the dot. The Landau levels in the accumulation layer 
become well-defined as the ratio T/lu c < 1, and their 
formation is observable in the current flowing through 
the quantum dot ground state by the presence of several 
peaks. The Landau levels are shifted towards the emit- 
ter Fermi level by the magnetic field. This is reflected in 
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fig. |3jby the fact that the peak associated to each Landau 
level shifts towards the threshold voltage as the magnetic 
field increases. 

As the overlap integrals increase with B, the oscilla- 
tions shown by the current as a function of B increase 
in amplitude with the magnetic field. It is important to 
notice that for a 2D emitter the increase in the amplitude 
of the current with B is much more pronounced than the 
one seen for 3D contacts. This is due to the fact that in 
3D the increase in the overlap between the wavefunctions 
in the lateral directions is counteracted by a decrease in 
the electron momentum along the direction of the cur- 
rent, as shown in the previous section. The latter effect 
does not occur in a 2D emitter and the change in the am- 
plitude of the current resonances as a function of the field 
is only due to the change in the overlap integral. This is 
a key difference between tunnelling from 2D and 3D con- 
tacts. The contribution to the current associated with 
each Landau level is directly proportional to the overlap 
integral between the Landau level and the quantum dot 
wavefunction. As the overlap with the first Landau level 
is the only one which tends to one when B increases, the 
contribution of the first Landau level is dominant. 

To the best of our knowledge, such an increase has 
never been observed and only 3D-0D and 1D-0D trans- 
port has been reported. Nevertheless our theory clearly 
predict this effect in quantitative terms and highlights the 
profound difference in the behaviour of the fine features 
associated to the quantum dot states as a function of the 
magnetic field depending on whether an accumulation 
layer is formed or otherwise. Furthermore, at least for 
the quantum dot ground state, our prediction is not af- 
fected by the presence of the electron-electron interaction 
in the quantum dot as when the ground state becomes 
resonant there is only one electron in the quantum dot. 

Finally, in fig. ^ we show the current threshold for dif- 
ferent dot states as a function of the magnetic field. As 
in the 3D-0D case (see fig. |J) we clearly observe here the 
selection rule over the angular momentum which sup- 
presses some peaks for the = 0,md > states, e.g.: 
peaks at about 23T for the current contribution of the 
first Landau level through the dot states irid = 1 and 2 
or peak at about 7T for the current contribution of the 
second Landau level through the dot state rrid = 2. Let 
us point out that the current contribution of the first 
Landau level at 23T through the n<j = 1, nid = is very 
weak but not suppressed as expected from the angular 
momentum selection rule. Furthermore the increase in 
the current peak with magnetic field can be observed for 
all those dot states with non-positive angular momentum 
and n = 0. 



V. CONCLUSION 

We have presented a theoretical model that allows us 
to calculate analytically the current flowing through a 
quantum dot state as a function of the magnetic field and 
the voltage bias. Calculations have been performed for 
3D-0D, 2D-0D tunneling. The results show that the am- 
plitude of the current is very sensitive to the angular mo- 
mentum of the quantum dot state. In particular the pres- 
ence of the magnetic field hinders tunnelling through the 
states with angular momentum co-linear with the mag- 
netic field. Furthermore in the case of tunnelling from 
2D accumulation layers, the current amplitude associated 
with certain dot states is predicted to increase hugely as 
an external magnetic field is applied. The analysis pre- 
sented in this work can be used to interpret transport 
experiments in quantum dots and directly access details 
of the wavefunctions of the zero-dimensional states Q . 
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FIG. 1. I-V for several dot states with a 3D emitter at T = 4K; V=0 corresponds to the threshold voltage. From the bottom 
to the top curve the magnetic field varies from to 10.25T with a step of 0.25T. The curves are shifted for clarity. 
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FIG. 2. Current maximum as a function of B for several dot states in the 3D-0D tunnelling case. 
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FIG. 3. I-V curves of the ground state for a 2D emitter at T = 4K; V=0 corresponds to the threshold voltage. From the 
bottom to the top curve the magnetic field varies from to 12T with a step of 0.5T.The curves are shifted for clarity. 
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FIG. 4. Current threshold as a function of B for different dot states in the 2D-0D tunnelling case. 
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